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ABSTRACT 

We examine r-process nucleosynthesis in the neutrino-driven wind from the thick accretion disk 
(or "torus" ) around a black hole. Such systems are expected as remnants of binary neutron star or 
neutron star - black hole mergers. We consider a simplified, analytic, time-dependent evolution model 
of a 3Mq central black hole surrounded by a neutrino emitting accretion torus with 90 km radius, 
which serves as basis for computing spherically symmetric neutrino-driven wind solutions. We find 
that ejecta with modest entropies (~ 30 per nucleon in units of the Boltzmann constant) and moderate 
expansion timescales (~ 100 ms) dominate in the mass outflow. The mass-integrated nucleosynthetic 
abundances are in good agreement with the solar system r-process abundance distribution if a minimal 
value of the electron fraction at the charged-particle freezeout, Y e ^ m i n ~ 0.2, is achieved. In the case of 
^,min ~ 0.3, the production of r-elements beyond A ~ 130 does not reach to the third peak but could 
still be important for an explanation of the abundance signatures in r-process deficient stars in the 
early Galaxy. The total mass of the ejected r-process nuclei is estimated to be ~ 1 x 1O~ 3 M0. If our 
model was representative, this demands a Galactic event rate of ~ 2 x 10~ 4 yr _1 for black- hole-torus 
winds from merger remnants to be the dominant source of the r-process elements. Our result thus 
suggests that black-hole-torus winds from compact binary mergers have the potential to be a major 
but probably not the dominant production site of r-process elements. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — binaries: close — stars: abun- 
dances — stars: neutron 



1. INTRODUCTION 

In the past decades, core-collapse supernovae (CC- 
SNe) have been considered to be the most promising 
astrophysical site for providing physical conditions 
suita ble for the rfrapid neutr on capture)-process (see, 
e.g., iThielemann et al.l 120111 for a recent review ) . 
The sce narios include the n e utrino-driven w i nd of 
CCSNe dWooslev et all Il99l iTakahashi et all Il99l 
Qian k Wooslevlll996t lOtsuki et alj|2000b IWanaio et aj 
20011: iThompson et al.l 120011 : iFarouai et all 1201 " 
prompt explosions of CCSNe (Sumiyo shi et al.l 120011 
or of electron-capture supernovae (ECSNe; a subset of 
CCSN e arising from collapsing oxygen-neon-magn esium 
cores. iHillebrandt et al.l 119841 : IWanaio et all 120031 ). and 
the shoc ked surface layers of the stellar core in collapsing 
ECSNe (INing et al.H2"00l . 

However recent hydrodynamical simulations of CC- i T^^l^^^T rr^F 
SNe (e.e.. iBuras et all 1200a IMarek k Jankall2009l) and 
of ECSNe pOtaura et al.l 120061 Uanka et al.l I2008D do 
not support the prompt explosion or the shocked sur- 
face layer scenarios. The nucleosynthesis calculations 
with one-dimensional hydrodynamical results of EC- 
SNe do not confir m the production of elements be- 
yond A ~ 90, either (jHoffman et alj|2008t IWanaio et all 
|2009| ). Furthermore, recent long-term simulations of 
CCSNe and ECSNe show that the neutrino-driven out- 
flows are proton-rich all the way (|Fischer et al.l I20T0I 
iHMepohl etHI f20Toh . which poses a severe difficulty 
to the scenario that has been favored for a long time, 
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neutrino-driven winds from proto -neutron stars (PNSs) 
in CCSNe and ECSNe. Recently, IWanaio et ail (|2011ali 
suggested on the basis of two-dimensional models of EC- 
SNe that a weak r-process could lead to the production of 
trans-iron elements in the early neutron-rich convective 
blobs of such SNe, but no heavier than A ~ 120. 

In contrast, another popular scenario of the as- 
trophysical r-proc ess, the mergers of double neutron 
stars (NS-NS, e.g..lRosswog et al.lll999HRuffert k Jaiikal 
19991: IShibata k UrvfJ 12000: Shi bata k Taniguchil 120061 
Oechslin et al.l 1200 Tt i or of black hole - neutron star bi- 
naries (BH-NS, e.g.. Uanka et al .lll999l: IShib ata k UrvfJ 
2006i IRuffert k Jankal 12010ft . has not been satisfac- 
torily explored. The decompression of dynamically 
ejected neutron-rich crust matter from NS-NS (or BH- 
NS) mergers was suggested to be an alternative or addi 
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, Meyer 1989; iFreiburghaus et al 
20051 1201 ID Hyper-massive NSs 
resulting immediately after NS-NS merg- 
Sekiguchi et"all l20lll iRezzolla et ail [201lt 
1201 If ) . giving rise to magneti- 
outflows for ~ 10- 
r-processed mate- 
both NS-NS 



(HMNSs) 

iBauswein k Jankal 
cally driven and neutrino-driven 
100 ms, are also suggest ed to eject 
rial ()Dessart et al.l I2009D . Furthermore, 



(after a possibly only short HMNS phase) and BH- 
NS (without a HMNS phase) mergers are expected to 
form a neutrino radiating accretion torus around the 
relic black hole, giving rise to neutrino-driven winds 
(hereafter, BH-torus winds) and potential short-duration 
gamma-ray-burst (GRB) jets. BH-torus winds are also 
expected t o provide suitable phys i cal conditions for the 
r-process (IRuffert k Jankal Il999l: ISurman et all 2008; 
iMetzger et al.ll2008bf ). ~ 
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Merger scenarios have been disfavored compared to 
those of CCSNe partly due to discrepancies between 
Galactic chemical evolution models and the spectro- 
scopic analyses of Galactic halo stars. The low Galac- 
tic ev ent rate (7 x 10 _6 -3 x 10~ 4 vr" 1 . IBelczvnski et al.l 
l2002f ) of mergers and the lon g lifetimes of binary sys tems 
(~ 1 Myr or 100-1000 Mvr. iBelczvnski et alll2002l) are 
expected to lead to the delayed appearance of r-elements 
in the Galactic history wi th too large star-to-star scatter - 
ing of their abundances (|QiarJl2000t lArgast et al.ll2004D . 
These facts seem to be in conflict with the observa- 
tional results of Gal ac tic halo stars (IHonda et al.l I2004t 
iFrancois et al.l I2007D . iBaneriee et al.l (|2011[ ) suggested 
that the early enrichment of r-elements might be due 
to neutrino-induced r-processing in the compact helium- 
shells of CCSNe of low-metallicity stars, and the con- 
tribution from mergers could follow only at a higher 
metallicity. However, recent studies of Galactic chemi- 
cal ev olution based o n the hierarchical clustering of sub- 
halos (jPrantzos 200(| also Y. Ishimaru 2011, in prepara- 
tio n) or with various binary popula tion synthesis mod- 
els (|De Ponder fe Vanbevere n 2004) do not exclude NS- 
NS and BH-NS mergers as the dominant astrophysi- 
cal site of the r-process in the early Galaxy. The rea- 
son of observed star-to-star scattering only found in r- 
elements (but not in a and ir on-group elements) has not 
been fully understood , either (llshimaru fe Wanaiolll999l: 
Tsuiimoto et alJ Il999t lArgast et al. | l2000t lArnone et al.1 
20051 iKarlsson fe Gustafssonl 120051: ICescuttil 120081 ). For 
these reasons, NS-NS and BH-NS mergers cannot be ex- 
cluded as the primary source of r-elements in the Galaxy. 
More studies of nucleosynthesis are highly desired. 

In this paper, we examine r-process nucleosynthe- 
sis in BH-torus winds, which are expected to be com- 
mon to both NS-NS and BH-NS mergers. There exist 
few previous studies of nuc l eosyn t hesis relevant to these 
conditions dSurman et al.1 120081: iMetzger etaLI 120091: 
ICaballero et al.ll2011l) . whi ch are based on param etrized 
outflow conditions. While IMetzger et al.l (120091) di scuss 
viscously driven mass eiecta. iSurman et al.1 (|2008T ) and 
ICaballero et al.l (|2011l ) examined some phenomenologi- 
cally chosen trajectories leading to at least a weak r- 
process, but did neither discuss time-dependences nor 
did provide the integrated abundance distribution or the 
ejected amount of r-process nuclei. Currently, multi- 
dimensional simulations of the wind phase after the for- 
mation of a stable accretion torus are not yet avail- 
able. Hence, we adapt a semi-analytic, spherically sym- 
metric, general relativistic steady-state wind model for 
nucleosynthesis calculations of BH-torus winds. The 
model has originally been developed for the studies 
of the r-process in neutrino-driven winds of CCSNe 
(|Wanaio et al.ll200TI ). 

Our paper is organized as follows. In § 2, we describe 
our model of stationary BH-torus winds and discuss basic 
outcomes derived from the wind solutions. In § 3, a phe- 
nomenological, time-evolutionary model of the neutrino 
luminosities of the torus is introduced, which is needed 
to determine the initial composition for nucleosynthesis 
and to calculate the mass-integrated yields as well as the 
ejecta mass. The results of nucleosynthesis calculations 
with the wind solutions are then presented in § 4, along 
with the mass-integrated abundances ejected from the 
BH torus. In § 5, we discuss the potential role of BH- 




FlG. 1. — Sketch of our model design for BH-torus winds. A 
rotating BH with the mass Mbh = 3Mq is located in the cen- 
ter of an accretion torus ("neutrino surface") that lies between 
2i?s and 10-Rg from the center, where ijg is the Schwarzschild 
radius (= 8.86 km). The wind is assumed to be radial. Thus 
we replace the complex neutrino-surface geometry by a spherical 
outflow model considering a neutrinosphere with radius i?„ = D 
around the gravity center (e.g., the asterisks on the dotted circle). 

torus winds as the origin of r-elements in the Galaxy. A 
summary of the paper follows in § 6. 

2. MODELING BH-TORUS WINDS 

Our treatment of BH-torus winds is based on the 
semi-analytic, spherically symmetric, general relativistic 
model of PNS winds that has bee n developed for stud - 
ies of r-process nucleosynthesis bv lWanaio et al.l (|2001l) . 
The mass of the central BH is taken to be Mbh = 
3M , which can be considered as representative of NS- 
NS mergers (or BH-N S mergers with a low-mass BH, 
Bclc zvnski et all 12008). As illustrated in Fi gure 1, the 
rotating accretion torus around a spinning BH is de- 
fined as the "neutrino surface" and is assumed to lie 
between 2R S (= 17.7 km) and 10i? s (= 88.6 km) from 
the center (where i?s = 8.86 km is the Schwarzschild 
radius) as suggested by detailed hydrodyn amical simu- 
lations of NS-NS and BH-NS merging (e.g., Uanka et all 
119991 : iRuffert fc Jankalll999l : iSetiawan et"alll2006f FI. 

In order to connect the asphcrical configuration of the 
BH-torus system to our spherical model for the wind out- 
flows, an arbitrary point on the torus is replaced by a 
point on a hypothetical neutrinosphere with the same 
distance from the center, R v = D (dotted circle in Fig- 
ure 1). The wind trajectory reaching away from the neu- 
trinosphere, which yields a description of the dynamical 
and thermodynamical outflow properties, is then derived 
in the same manner as for spherical PNS winds by solv- 
ing the general relativistic stationary equations of mass, 
momentum, and e nergy conservation (Eqs. (l)-(3) in 
IWanaio et al.|[200ll ). The Schwarzschild geometry due to 
the central BH is included here, but the gravitational ef- 

3 The radius of the inner tip (or innermost stable circular orbit) 
is smaller for a BH co-rotating with the torus than for the non- 
rotating case. Our case, the innermost radius of D = 2R B , corre- 
sponds to a BH spin parameter of ~ 0.6 (Fig.l in Sctiawan et al. 
2006). The BH rotation is expected not to lead to major differ- 
ences in the dens ity or temperature d istributions except for the 
innermost region (Sctiawan et al. 2006 ) but to a moderate increase 
of the torus mass (~ 50%, Sctiawan et al. 2006). Because the in- 
nermost winds make an only small contribution to the integrated 
abundances, and we consider the case with M e j <C A/torus (cf. , 
Eqs. (11) and (16)), the effects of BH rotation are expected to be 
secondary. 
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Fig. 2. — Assumed neutrino luminosities (L v ) as functions of 
the distance from the center (D). L v is assumed to increase lin- 
early from L„,o/100 to L u fi (denoted in the legend in units of 
10 51 erg s" 1 ) between 2R S (= 17.7 km) and 5fts (= 44.3 km; 
dotted line) and adopt a constant value at larger distances D. 

fects of the torus and wind masses are very small and thus 
neglected. Rotation of the mass-losing object is neglected 
as well. The equation of state for ions (ideal gas) and ar- 
bitrarily degenerate, arb itrarily relativistic e lectro ns and 
positrons is taken from iTimmes fc Swestvl (|2000). The 
average neutrino energies are taken to be 15, 20, and 
30 MeV, for electron neutrino, electro n antineutrino, and 
heavy- lepton neutrin os, respectively (jJanka et alj|1999t 
iSetiawan et al.l l2006[ ) . The neutrino luminosities of all 
the flavors are assumed to have the same value L v in our 
wind model. The electron fraction (number of protons 
per nucleon), Y c , is assumed to be 0.5 (for a justification, 
see0). At the inner boundary, the density is taken to be 
p = 10 10 g cm~ 3 (a different choice has no big impac10), 
and the temperature, T, is taken such that neutrino heat- 
ing and cooling balance each other (~ a few MeV). The 

velocity (or equivalently, mass ejection rate, M) at the 
neutrinosphere is determined such that the wind becomes 
supersonic through a sonic point. 

As anticipated from Figure 1, the neutrino flux from 
the outer regions of the torus is shielded in the vicin- 
ity of the BH by the presence of the torus itself. In 
order to mimic this effect in our spherical models, we 
simply assume that L v increases linearly from L^o/100 
to L vfi between 2R S (= 17.7 km) and 5R S (= 44.3 km) 
and adopts a constant value on the outer side. Figure 2 
shows the assumed profiles for selected L^o in units of 
10 51 erg s _1 . We define the outflows from R v < 5Rs 
and R v > 5Rs as the "inner" and "outer" winds, respec- 
tively. The calculated radial profiles of velocity it, density 
p, and temperature T for the cases with R v = 2i?s (solid 

4 As described in § 3, L Vll T <C L„ e < La a and Y c <C 0.5 are 
expected in the early phase of BH-torus outflows. However, as- 
suming L v = ^(ii/ c + Lv c ) for all neutrino kinds and Y B = 0.5 
is reasonably good for computing the dynamics of BH-torus winds 
because the wind-driving energy deposition by v B and v e is very 
similar while that of heavy-lepton i/s is small. Note that these 
assumptions for L v and Y e are applied only in this section. 

5 The wind profile hardly depends on the boundary density for 
> 10 10 g cm"'. Tests show that the wind profiles remain essen- 
tially unchanged when the inner boundary density is varied be- 
tween values somewhat above 10 9 g cm -3 (where the density gra- 
dient of the wind changes; Fig. 3, left-middle panel) and below the 
neutrinospheric value (around 10 11 g cm -3 or higher). 



lines), 5Rs (dashed lines), and 10i?s (long-dashed lines) 
are shown in Fig. 3 (left panels) in dependence on the dis- 
tance r from the center. L„.o is taken to be 10 53 erg s _1 . 
Right panels in Fig. 3 display the temporal evolution 
of r, p, and T of selected ejected mass elements, where 
the time coordinate t — is set at the inner boundary 
(p = 10 10 g cm~ 3 ). 

The above assumption of the neutrino luminosity pro- 
files (Fig. 2) implies fundamental differences in the neu- 
trino heating properties between the inner and outer 
winds. In Figure 4, the neutrino heating and cooling 
rates, q (in units of erg g _1 s -1 ), for the cases with 
R v = 2Rs (solid lines), 5i?s (dotted lines), and 10i?s 
(dashed lines) are displayed as functions of the radial dis- 
tance from the center, r. L^o is taken to be 10 53 erg s . 
Heating is due to v e and v e captures on free nucleons 
(vN), neutrino scattering by e~ and e + {ve), and vv 
pair annihilation to e~e + pairs (yv). Cooling is caused 
by e~ and e + captures on free nucleons (eN) and e~e + 
pair annihilation to vv pa irs (ee). All the rates are taken 
from lOtsuki et al.l (|2000l . see their Eqs.(8)-(16)), where 
the gravitational redshift of neutrino energies and the 
bending of trajectories due to general relativistic effects 
are fully taken into account. 

As can be seen in Figure 4, q V N dominates the heat- 
ing rates in the outer winds (dotted and dashed lines). 
The rate q ue plays only a secondary role, and q vv is 
negligible. This can also be foun d in the studies of 
PNS winds (IQian fc Wo oslevl [l99l l Otsuki et all llM 
iThompson et al.l 120011 ). In contrast, in the innermost 
winds (solid lines), q vv plays the dominant role for heat- 
ing. This is a consequence of the substantially smaller 
mass ejection rate, M (Fig. 5; defined as that from the 
corresponding neutrinosphere), driven by the small L v 
from the innermost region of the torus (Fig. 2). This 
leads to very small va lues of p and als o of T in the vicin- 
ity of the torus (cf. iWanaiol l2006bl . for a PNS case). 
For a given radius r, the heating rates scale with L v , 
p, and T according to q„N oc L v , q ue oc L y p _1 T 4 , and 
q vv oc L v 2 p~ x . As a result, the reduction of p and T at 
small r boosts q vv due to radially increasing L v (as given 
in Fig. 2) to much higher values than the other r ates (see 
a simi lar discussion for anisotropic PNS winds, IWanaiol 
2006b). Indeed, this has been discussed as a promising 
mechanism to power short GRB j ets presumably arising 
from NS-NS or B H-NS merging (|Ruffert fc Jankalll99a 
IJanka et al.| [l999. and references therein). 

The greater q vv in the inner winds leads to substan- 
tially faster (terminal) radial velocities, u, higher asymp- 
totic entropies per nucleon, S, and short expansion 
timescales, r (defined as the e-folding time of T below 
0.5 MeV). As shown in Figure 6, the innermost winds be- 
come relativistic, achieving u ~ (0.3-1) c (where c is the 
velocity of light). Figures 7 and 8 also show that the in- 
nermost winds obtain extremely high S (~ 1000 in units 
of the Boltzmann constant, &b) and short t (~ 1 ms). 
In contrast, the outer winds have modest u (~ 0.1c), S 
(~ a few 10 fes) and r (~ a few 100 ms). For a given D, 
a larger L v p leads to higher u and shorter r because of 
the greater neutrino energy deposition. However, the re- 
sulting larger M reduces the neutrino heating per mass, 
leading to slightly smaller S. 

Note that the three-dimensional simulations of BH ac- 
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Fig. 3. — Left: Radial profiles of velocity (top), density (middle), and temperature (bottom) for the cases with ij„ = 2i?g (solid lines), 
5i?s (dashed lines), and 10i?s (long-dashed lines) versus distance r from the center. L v fi is taken to be 10 53 erg s _1 . Asterisks mark the 
sonic points. Right: Temporal evolution of radius (top), density (middle), and temperature (bottom) of ejected mass elements for the same 
cases. Time is set to t = at the inner boundary (p = 10 10 g cm -3 ). 

the inner boundary to the heating region where q peaks 
(Fig. 4). Once the matter reaches this maximal- heating 
region, where the density is already significantly lower, 
the subsequent matter motion would not be substantially 
disturbed by the time-evolving torus conditions. We thus 
regard our steady-state model applicable if T; n j <C 100 ms 
is satisfied. Figure 9 shows the injection distances, de- 
fined as di n j = r(q = q max ) — D, as functions of D for se- 
lected L v ^. The inner (except for innermost) winds have 
relatively small c?j n j (~ 10 km), while the outer winds 
have greater di n j increasing with D (from 10 km to 30- 
40 km). Note that di n j is relatively inse nsitive to L^.q as 
also found in the PNS case (see Fig. 4 in I Thompson et aLI 
120011 ). Figure 10 displays the injection timescales Tj n j 



cretion tori by Sctiawa n et al.l (|2006| ) indicate that the 
matter in the inner region moves quickly towards the 
central BH with short accretion timescales of ~ 1-10 ms 
(with the radial velocities of ~ (0.01-0. l)c). Moreover, 
the lifetime of an accretion torus is estimated to be of 
the order of ~ 100 ms. This is substantially shorter 
than the period of a PNS wind phase (~ 10 s). It 
is important, therefore, to test the applicability of our 
steady-state wind solutions to the BH-torus case. More 
specifically, the matter should escape from the BH-torus 
surface region within a sufficiently short period of time 
(< 100 ms). 

To test this, we introduce the "injection timescale" 
Ti n j for each wind, during which the matter moves from 
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Fig. 4. — Neutrino heating and cooling rates, q, for the cases 
with i?„ = 2i?g (solid lines), 5-Rg (dotted lines), and 10i?s (dashed 
lines) as functions of radius from the center, r. L„ o is taken to 
be 10 53 erg s — 1 . Heating is due to v e and v £ captures on free 
nuclcons (vN), neutrino scattering by e~ and e+ {ve), and vv pair 
annihilation to e~ e+ pairs (vv). Cooling is due to e~ and e+ 
captures on free nucleons (eN) and e~ e+ pair annihilation to vv 
pairs (ee). The net total rates are indicated in black. The vertical 
dashed line indicates D = 5-f?s- Note that L v increases with r as 
in Fig. 2 for the i?„ = 2i?,g case and that the steep drop of q„jv 
for this case is due to the decreasing free nucleon abundance by 
a-particle formation. 
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Fig. 5. — Mass ejection rates, M, from the equivalent neu- 
trinosphere R v = D, as a function of D for L„,o (in units of 
10 51 erg s — 1 ) denoted in the legend. The vertical dashed line 
indicates D = 5i?s- 
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Fig. 6. — Terminal radial velocities as functions of D for L v fi (in 
units of 10 51 erg s — 1 ) denoted in the legend. The upper bound is 
the velocity of light. The vertical dashed line indicates D = 5R$. 
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Fig. 7. — Asymptotic entropies as functions of D for L„ o ( m 
units of 10 51 erg s" 1 ) denoted in the legend. The vertical dashed 
line indicates D = 5i?g. 
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Fig. 8. — Expansion timescales (the e-folding times of temper- 
ature below 0.5 MeV) as functions of D for L v fi (in units of 
10 51 erg s — 1 ) denoted in the legend. The vertical dashed line 
indicates D = 5Rs. 
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Fig. 9. — Injection distances, di n j, as functions of D for L„,o (in 
units of 10 51 erg s _1 ) denoted in the legend. Each wind reaches 



the radius with the maximal q at r 
line indicates D = 5-Rg . 



D + dj n j . The vertical dashed 



(solid lines) as functions of D for selected L^ q. We find 
particularly large r; n j for the innermost winds (see also 
right panels in Fig. 3). It should be noted, however, that 
the injection timescales could be overestimated owing to 
the ambiguity of the boundary density (see footnote 5) 
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Fig. 10. — Injection timescales, Tj n j, as functions of D for L„,o 
(in units of 10 erg s — 1 ) denoted in the legend. Each wind reaches 
the radius with the maximal q at t = r; n j, where t is set to at 
r = D (solid lines) or at r = D + 1 km (dotted lines; see text). 
Dashed lines indicate the elapsed times from r = D + 1 km to the 
sonic points. The vertical dashed line indicates D = 5i?s . 

and thus of the boundary velocity (cf., M = A-nr 2 pu). 
To avoid this, the elapsed times from r = D + 1 km to 
r = D + dj n j are also drawn in Fig. 10 (dotted lines). 
This indicates that the large r- m j for the innermost winds 
is largely due to the small initial velocities (see Fig. 3, 
top-left panel) arising from our choice of the boundary 
density (p = 10 10 g cm~ 3 for all the cases). 

Considering the condition Tj n j <C 100 ms, we find that 
our steady-state treatment can be (at least marginally) 
applicable for the winds with L^o > 50 (in units of 
10 51 erg s _1 ). Note that, once the matter reaches 
r = r(g max ), the outflow becomes supersonic within a 
short period of time. The supersonic outflows are soni- 
cally disconnected from the surface and are not affected 
by the time-evolving BH-torus any more. Dashed lines in 
Fig. 10 indicate the elapsed times from r = D + 1 km to 
the sonic points, which are ~ 2r; n j except for the inner- 
most region. The total ejecta mass might be somewhat 
overestimated owing to the longer r in j for lower L„, - As 
we find in Fig. 5, however, M for L v $ < 50 (in units of 
10 51 erg s _1 ) is substantially smaller and its contribution 
to the total mass will be sub-dominant. The ejecta mass 
from the inner region might be also overestimated be- 
cause of the fast accretion timesca les estimated from hy - 
drodynamical studies (~ 1-10 ms. lSetiawan et aLl l2006). 
As we will see below, however, the ejecta from the outer 
region dominate the total amount (~ 90%; § 3). It is 
also important to note, taking our injection timescale 
condition seriously, that the outer-torus region beyond 
^100 km is not expected to contribute significantly to 
the total ejecta mass because of long injec tion timescales. 

It s hould be noted that the work by ISetiawan et al.l 
(2006) was based on a neutrino trapping/leakage scheme 
and employed a simple a-model for the gas viscosity. The 
dynamics of outflows from BH tori will be highly depen- 
dent on neutrino transport as well as magnetic fields, 
which are to be explored in f uture studies of hyper ac- 
creting BHs. This did not allow lSetiawan etaTI (|2006D to 
describe mass loss from the tori by energy transfer and 
deposition associated with neutrino transport, nor did 
their models adequately account for magnetically-driven 
mass loss. Outflows and winds were previously discussed 



bv iPruet et all (12001: iMetzger et all |2008allbl. [2 009); 
iSurman et aTT(|2008l ) : TCaballero et al.l ()201lD : iRosswod 
(|2011f ). although all of those works did not employ self- 
consistent hydrodynamical models of this phenomenon. 

3. TIME EVOLUTION OF THE BH-TORUS 

3.1. Why is the BH-torus wind neutron-rich? 

As mentioned in § 1, recent hydrodynamical simula- 
tions of_CCSNepmh£t th^ is proton- 
rich ()Fischer et a lj[2010l:IHu depohl et al.ll2010h. The BH- 
torus wind from "collapsars" (MacFad ven fc Wooslevl 
1999) resulting from collapsing rapidly rotating stel- 
lar cores is also suggested to be proton-rich (e.g., 
IKizivat et all [20101) . The BH-torus from NS-NS and 
BH-NS mergers, originating from decompressed, very 
neutron-rich NS matter, is different from the collapsar 
case. The initially neutron-rich torus leads to more pro- 
duction of v e by e + captures on neutrons, 

n + e + — > v c + p, (1) 

than v e production by e~ captures on protons, 

p + e~ — > v c + n, (2) 

and thus one obtains L Pc > L Ve . This allows the BH- 
torus wind (outside of the torus at lower densities) to 
remain neutron-rich owing to more D e captures on free 
protons, 

D c + p — > n + e + , (3) 

than v e captures on free neutrons, 

v c + n — > p + e~ . (4) 

However, despite the fact that the neutrinospheric lay- 
ers in a PNS are also neutron-rich, the simulations cited 
above show that PNSs produce essentially equal lumi- 
nosities of v c and D e , resulting in the proton-richness 
of the ejecta because of the neutron-proton mass differ- 
ence, which reduces v c captures compared to v c captures. 
So what exactly is the difference between the PNS and 
merger BH-torus cases? This can be explained as follows. 

For PNSs, the main point is the fact that the PNS is in 
a phase of neutronization. It evolves from more proton- 
rich (symmetric) initial conditions to more neutron-rich 
final conditions. In course of this transition, its core 
radiates more v c than v c (see Eqs. (1) and (2)). On 
their radial way out of the PNS through highly neutrino- 
opaque layers this initial v c flux is converted into equal 
u e and u c fluxes by multiple absorption and reemission 
processes (and p and r fluxes of similar size). The es- 
caping v c and v e luminosities are roughly equal because 
the neutrinospheric layers with their typical densities 
(~ 10 11 g cm -3 initially and up to ~ 10 13 g cm -3 at very 
late times) reach their so-called "/3-equilibrium" state 
(defined by equal number production of v c and u e and 
Y c = const, in time) within some tens of milliseconds up 
to about ~ 100 ms. This is much shorter than the PNS 
cooling timescale of seconds. On the long PNS evolution 
timescale the emission is therefore mostly characterized 
by essentially equal L„ c and Lp a (with a small difference 
making sure that there is a net i/ c -flux out of the still 
neutronizing high-density PNS core). 

In contrast, the hot post-merger torus is composed of 
decompressed neutron star matter. It starts out from a 
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very neutron rich initial state and (on its way to a new (3- 
equilibrium condition) protonizes gradually (Eq. (1)) be- 
cause of overall lower densities than it had initially. Since 
the density and temperature of the torus (on average 
some 10 11 -10 12 g cm -3 and about 3-10 MeV) are very 
similar to the neutrinospheric conditions of a PNS, the 
protonization of the torus proceeds on the same timescale 
as the neutronization of the neutrinospheric region in the 
PNS case. As mentioned above, this typical timescale is 
tens of milliseconds at the given densities, presumably 
up to ~ 100 ms. This is exactly a major fraction of the 
accretion timescale of most of the torus mass into the 
BH (and is the timescale for the protonization assumed 
in our simple model described in § 3.2). Because of the 
protonization, L„ e is higher than L„ e , in inverse analogy 
to the fact that L„ c is higher in the case of a collapsing 
stellar core during the early times when the postshock 
layer is still on its way to reach the /3-equilibrium condi- 
tion, starting out from a nearly symmetric initial state. 
For this reason, the neutrino-driven ejecta from the BH- 
torus are expected to be neutron-rich (see Eqs. (3) and 
(5)) for a significant duration of the wind phase. 

Note that a direct consequence of the fact that the 
torus conditions in density and temperature are roughly 
similar to the neutrinospheric conditions in a PNS is the 
finding that the /i and r neutrino luminosities are 5-10 
times lower th an the v P . and v P luminosities in the cas e 
of the tori (see lJanka et al.l ll999: Ru ffert fc Ja nka 1999). 
In a PNS the /x and r neutrinos are produced at much 
higher densities than u e and v e . Such densities are not 
reached by the accretion tori. 

3.2. Time Evolution of Neutrino Luminosities 

In order to calculate the nucleosynthesis for each wind 
trajectory, the Y c value of the torus, which defines the 
initial composition, should be specified. In addition, the 
ejecta mass for each wind trajectory should be deter- 
mined to calculate the mass-integrated nucleosynthctic 
abundances. For this purpose, we assume that the time 
evolution of L v and L p can be written as 



where L,, 



L v< , (t) — L^i 



L Pe (t)=L v ,i 

,53 



-0 



1 + 2[ u 



(5) 
(6) 



10 erg s is chosen as a representative 



value of the initial neutrino luminosity at t = £i = 10 ms 
(e.g. . Uanka et al. 1999; Scti awan et al.ll2006[) . The power 
exponent p is taken to be 1.3 (jMetzger et al.l I2008bt 
iLee et ahl 120091 ) . The value of a is uncertain and for 
simplicity assumed to be a modest number, a = 1, be- 
cause Lv c > L Ue is expected during a significant time 
of the torus evolution (§ 3.1). These relations give 



where iVb and (e„) are the total baryon number of 
the torus and the mean neutrino ener gy (~ 18 MeV, 
Uanka et al.l 119991 : iSetiawan et al.l I2006D averaged for v c 
and P c . This can be integrated to yield 



Y e (t)=Y e . 



Y P . 



2L V 



N h (e u ) 

2L V \ ti 



-(a+(9) 



dt' 



(a + 0-l)N h (e u ) 



-(a+/9-l) 



where Y c i = 0.05 is assumed for t he initial electron 
fraction (e.g., iRuffert &: Jan ka 1999). Further assum- 
ing that Y e asympto tes to Y e ,{ — 0.50 for t >• ti (e.g., 
iMetzger et"aT1(2008b[) . we get 



2L L/ i ti 



{a + p-l)N h {e v } 
This simplifies Eq. (8) to become 



Y e (t) = Y e ,f - (Y e , { - y 6)i ) ( - 



Y e ;. 



-("+/3-1) 



(9) 



(10) 



Eq. (10) provides Y c of the torus at a given evolutionary 
time, t. Note that the value of iVb (not made use of in 
this study) can be estimated from Eq. (9) and leads to 
the approximate torus mass, 



= N h m u 



2L y i ti m u 



(a + /3-l)(Y e , f -y e , i )(e 1/ > 



Q.1NM,) 



which is consistent with hydrody namical res u lts of 
compact binary me rgers (e.g., Uanka et al.l 119991 : 
IRuffert fc JMika][T99a . 

3.3. Connecting Steady Wind and BH-Torus Evolution 

Models 

Our description of BH-torus winds in § 2 is based on 
steady-state solutions for constant L„.o and Y c . In order 
to connect the computed wind trajectories to the time- 
evolving torus model described in § 3.2, we assume 



L*. (*)+£».(*) 



(12) 



This is a reasonable assumption because the wind from 
the outer torus is driven predominantly by v c and v e cap- 
tures on free nucleons (Fig. 4), whose heating effect is 
nearly symmetric between v c and v c . This is not true 
for the innermost torus, which however does not sig- 
nificantly contribute to the mass-integrated nucleosyn- 
thetic abundances as we will see later. Eq. (12) gives, 
with Eqs. (5), (6) and (10), L„ c , L Po , t, At, and Y e (as 
the initial value for nucleosynthesis) for a given L v $ in 
our steady-state model of BH-torus winds. These val- 
ues are listed in Table 1 for L„ o = 200, 190, • • • , 10 (in 
units of 10 51 erg s _1 ). Here, At are the time intervals 



L„ e /L Ve = 3 at t = tj, being in agreement with hy- At,- = tj+i - tj between the j-th and (j + l)-th wind 



drody namical results (jJanka et al.lll999t ISetiawan et al.l 
120061 ). and L Po = L„ e for t » h (jMetzger et al.ll2008bf) . 
Eqs. (5) and (6) yield an approximative time evolution of 
(the torus-averaged value of) Y c in the protonizing torus 
(§ 3.1; Eqs. (1) and (2)), 



trajectories. For the 20th wind trajectory (last line in 
Table 1), AHs determined such that the total evolution- 
ary time becomes t = 100 ms. 

The mass ejection rate for the j-th value of L v $ from 
the torus region Di < D < Di + i is calculated as 



Y. 



L p ,(t) - L v .(t) 
N b (e v ) 



2L,. 



t \ - {a+ ^ 



(7) 
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M, 



AttD 



2D? 
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TABLE 1 
Time-Evolution of the BH Torus 



j 








t [ms] 


At [ms] 




1 


200 


100 


300 


10.0 


0.290 


0.050 


2 


190 


96.4 


284 


10.3 


0.316 


0.066 


3 


180 


92.6 


267 


10.6 


0.346 


0.083 


4 


170 


88.9 


251 


11.0 


0.381 


0.100 


5 


160 


85.0 


235 


11.3 


0.422 


0.118 


6 


150 


81.0 


219 


11.8 


0.471 


0.135 


7 


110 


77.0 


203 


12.2 


0.530 


0.153 


8 


130 


72.9 


187 


12.8 


0.603 


0.172 


9 


120 


68.6 


171 


13.4 


0.692 


0.191 


10 


110 


64.3 


156 


14.1 


0.806 


0.211 


11 


100 


59.8 


110 


14.9 


0.953 


0.231 


12 


90 


55.1 


125 


15.8 


1.15 


0.252 


13 


80 


50.3 


110 


17.0 


1.42 


0.274 


14 


70 


45.3 


94.7 


18.4 


1.81 


0.296 


15 


60 


40.1 


79.9 


20.2 


2.41 


0.319 


16 


50 


34.7 


65.3 


22.6 


3.40 


0.344 


17 


40 


28.9 


51.1 


26.0 


5.26 


0.370 


18 


30 


22.7 


37.3 


31.3 


9.56 


0.398 


19 


20 


16.1 


23.9 


40.8 


24.7 


0.428 


20 


10 


8.68 


11.3 


65.6 


34.4 


0.461 


a in 


units of 10 51 


erg s — 1 . 











initial value for nucleosynthesis. 

where the torus is imagined to be a disk, ejecting matter 
perpendicularly toward the north and south directions. 
Here, Mj j is the mass ejection rate from the correspond- 
ing neutrinosphere with R v = Di (Fig. 5) for the j-th 

value of L„,o, and the mass flux density (M itj /4nDf) of 
the spherical wind calculation is weighted by the effec- 
tive torus surface element in Eq. (13). The mass ejection 
rate during time interval Atj from the entire torus with 
L Vt a = (L„fl)j is then given as 



= £Ar 



(14) 



Eq. (13) with the At from Table 1 gives the ejecta mass 

Am.ij — ArhijAtj (15) 

for the wind at D = Di and with L U fi — (L Ui o)j. From 
Eq. (15) the total ejecta mass in the BH-torus outflow 
during the first 100 ms is calculated as 



i.96 x w~ 6 m. 



0! 



(16) 



where the inner region (D < 5i?s) contributes only 
~ 10% (= 2.40 x 10~ 4 M Q ). The total ejecta mass in 
Eq. (16) is only a few percent of our assumed (initial) 
BH-torus mass (~ 0.1M©; Eq. (11)). This may justify 
our assumption of steady-state neutrino-driven outflows, 
provided that the lifetime of the BH-torus is > 100 ms, 
which can be expected to be the case, considering typical 
mass-accretion rates of BH-tori in hydrodyn amical mod- 
els (see, e.g., Fig. 4 in S etiawan et al.ll2006t ). Note that 
the initial high mass-loss phase (with M > 0.1M Q s -1 ; 
Fig. 5) for the outer winds lasts for only a short period 
of time (a few ms; Table 1). 

From Eq. (15) and with the numbers of Table 1, the 
ejecta mass Am e j (normalized by m C j) as a function of 
u can be constructed in a binned form as histogram 
(Fig. 11). We find that the sub-relativistic winds with 
u ~ 0.1c dominate the ejecta. The relativistic ejecta from 
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Fig. 11. — Histogram of the ejecta mass distribution Am c j (nor- 
malized by the total ejecta mass, m e j ) in linear (left vertical axis) 
and logarithmic (right vertical axis) scales as a function of the ter- 
minal expansion velocity u (in logarithmic scale). 
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Fig. 12. — Same as Figure 11, but as a function of the asymptotic 
value of the ejecta entropy S. 
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Fig. 13. — Same as Figure 11, but as a function of the expansion 
timescale r. 

the innermost torus are essentially unimportant. While 
the entropies are high, S ~ 100-1000 fee, in the inner- 
most winds (Fig. 7), the ejecta with modest entropies, 
S ~ 30fce, dominate (Fig. 12). Moreover, the winds with 
modest t ~ 100 ms dominate, and those with short r ~ 
1-10 ms make a very small contribution (Fig. 13). 

Note that L v ^ as well as (5 in Eqs. (5) and (6) would 
strongly depend on the mass of the torus as well a s 
the viscous parameter (Fig. 10 in Setia wan et al.l 2006). 
whose values are highly uncertain. A higher value of L v i 
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will increase the ejecta mass m e j in Eq. (16) as antici- 
pated from Fig. 5. A smaller (3 increases the ejecta mass 
as well (e.g., by about a factor of two when f3 = 0.7, 
roughly half of the original value, is chosen) because 
of a higher L v o for a longer time (leading to a higher 
mass ejection rate; Fig. 5). However, the nucleosyn- 
thetic abundance distributions would not be significantly 
changed. This is due to the fact that Eq. (10) can be 
written as Y c {t) — Y c j — (Y C: f — Y % i)L VB / L Vt i, i.e., as a 
function of L„ a only (for a — 1), being independent of 
/3. Because the change of T„ e is similar to that of L v p 
(Table 1), the Y c value at a given L^o is quite simi- 
lar, e.g., AY C < 0.02 for f3 — 0.7 is associated with the 
somewhat shallower decline of the neutrino luminosities. 
Accordingly, the integrated abundance distribution will 
be mostly independent of /3. 

4. NUCLEOSYNTHESIS IN BH-TORUS WINDS 

The nucleosynthetic yields in each wind trajectory are 
computed by solving an extensive nuclear reaction net- 
work code. The network consists of 6300 species between 
the proton and neutron drip lines, all the way from sin- 
gle neutrons an d protons up to the Z = 110 isotopes (for 
more detail, see iWanaio et al.ll2011bf ) . Neutrino interac- 
tions on free nucleons and a-particles are also taken into 
account, but fission reactions are not included. As will be 
described in § 4.3, the effect of (neutron-induced) fission 
during the r-process is expected to be of minor impor- 
tance in our explored cases. Our fiducial case ("case 1") 
is defined by initiating each nucleosynthesis calculation 
when the temperature decreases to Tg = 10 (where Tg is 
the temperature in units of 10 9 K), at which nuclear sta- 
tistical equilibrium (NSE) is immediately recovered from 
arbitrary initial compositions. The initial composition 
we adopt is Y c and 1 — Y c for free protons and free neu- 
trons, respectively, using the Y c value of the torus for a 
given L V fi (see Table 1). The calculations are carried out 
for wind solutions with L v ^ = 10-200 in steps of 10 (in 
units of 10 51 erg s _1 ; Table 1) and D = (2-10) R s = 
17.7-88.6 km with intervals of 0.1i? s = 0.886 km, i.e., 
for 1620 trajectories in total. 

4.1. Electron Fractions 

Figure 14 shows the variation of Y c and of the mass 
fraction of a-particles, X a , as functions of (decreasing) 
temperature for L v ,g = 2 x 10 53 erg s _1 and D = 5Rs 
(case 1; solid lines). This wind trajectory has the low- 
est initial (torus) Y c of 0.05 (indicated by a horizontal 
solid line) of our considered models. As can be seen in 
Fig. 14, Y c increases during the early wind expansion 
phase from its initial value (at Tg = 10; vertical solid 
line) to 5^,2.5 (at Tg = 2.5, considered as corresponding 
to the onset of the r-process phase; vertical dashed line). 
This is caused by several effects. The first effect is that, 
as soon as the nucleosynthesis calculation starts, Y c re- 
la xes toward the equilib rium value, Y c ^ a (e.g., Eq. (77) 
in iQian fc Woo slev 1996|), due to v e and v c captures on 
free nucleons (Eqs. (3) and (4)). The asymptotic value 
for this wind trajectory is Y c . a = 0.21. The second ef- 
fect is associated with the continuous a-particle forma- 
tion (peaking at Tg m 7 in Fig. 14) from recombination 
of free protons and free neutrons produced by neutrino 
capture reactions (E qs. (3) and (4)). Th is drives Y e to- 
ward 0.5 ("a-effect": lMever et atlll998al ). As a result of 




Fig. 14. — Evolution of the electron fraction (Y c ; red) and of the 
mass fraction of a-particles (X a ; blue) as functions of (decreasing) 
temperature for case 1 (solid lines; only for Tg < 10) and case 2 
(dotted lines). T„.o = 2 X 10 53 erg s — 1 and D = 5i?s are taken for 
both cases. For case 2, three cases of the Y e evolution are shown 
(for Tg > 10), in which all four reactions (Eqs. (1)— (4)), only i/ e 
and P e captures on free nucleons (Eqs. (3) and (4)), and only e~ 
and e+ captures on free nucleons (Eqs.(l) and (2)) are included 
(indicated by "all", V, and "e", respectively). The horizontal 
and vertical solid lines indicate, respectively, the initial (torus) Y e 
(= 0.05) and Tg = 10 (at which the computation is switched to 
making use of the full reaction network). The vertical dashed line 
at Tg = 2.5 indicates where r-processing begins. 



these combined effects, Y c 2 .5 (= 0.18) is sizably higher 
than the initial value of 0.05. 

It is important to note that v c and D c captures on free 
nucleons (Eqs. (3) and (4)) and also e~ and e + captures 
on free nucleons (Eqs. (1) and (2)) will also operate be- 
fore the temperature decreases to Tg — 10. We therefore 
consider another case ("case 2"), in which each nucle- 
osynthesis calculation is started from closer to the neu- 
trinosphere (at R y = D = 5i?s and Tg = 60 in this 
case) including the four reactions of Eqs. (l)-(4) (la- 
belled "all" in Fig. 14; dotted lines). The computation 
is then switched to making use of the full reaction net- 
work at Tg — 10. Fig. 14 also shows the results when 
including only v c and v c captures ( ' V ; Eqs. (3) and 
(4)) and e~ and e + captures ("e"; Eqs. (1) and (2)) for 
Tg > 10. We find that the Y e immediately jumps to 
~ 0.15 mostly by neutrino captures, continues to gradu- 
ally increase to 0.26 by competition of all four reactions 
of Eqs. (l)-(4), and decreases slightly toward Y e a (but 
still Y c = 0.23 > r c , a = 0.21 at Tg = 10). The increase of 
Y e for Tg < 10 is solely due to the a-effect. As a result, 
Y c = 0.29 at T 9 = 2.5, which is ~ 0.1 higher than that in 
our fiducial case 1. 

Figure 15 shows the Y e values when the temperature 
has decreased to Tg = 2.5, 1^,2.5, over the entire D range 
for selected L v $ cases (denoted in the legend). The shifts 
of Ye from the initial values (see Table l) are particularly 
prominent in the inner winds, where the matter experi- 
ences a strong a-rich freezeout owing to the high 5" and 
short t (Figs. 7 and 8; cf. Fig. 17). As described above, 
the Ye.2.5 values are significantly different between cases 1 
and 2, in particular for the high L^o (i.e., low Y e ) winds. 

We emphasize, however, that this is based on our 
simplified wind models and the exact Y c evolution for 
Tg > 10 would be highly dependent on the detailed den- 
sity structure as well as the highly uncertain neutrino 
field in the vicinity of the torus, where only transport 
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Fig. 15. — Electron fractions for case 1 (solid lines) and case 2 
(dotted lines) when the temperature has decreased to 2.5 X 10 9 K 
as functions of D for Z/„ o ( m units of 10 erg s — 1 ) denoted in the 
legend. The vertical dashed line indicates D = 5Rs- 
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Fig. 17. — Mass fractions of a-particles for case 1 (solid lines) 
and case 2 (dotted lines) as a function of D for L„,o (in units 
of 10 51 erg s _1 ) denoted in the legend. The vertical dashed line 
indicates D = 5R$ ■ 
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Fig. 16. — Neutron-to-seed ratios (measured at Tg = 2.5) for 
case 1 (solid lines) and case 2 (dotted lines) as functions of D for 
L V: o (in units of 10 51 erg s _1 ) denoted in the legend. The vertical 
dashed line indicates D = 5Rs- 

calculations can give an appropriate description of the 
energy and direction distribution of the neutrinos. We 
therefore consider the Y c 2 .5 differences between cases 1 
and 2 as possible uncertainties of Y c in this study. As 
can be seen in Fig. 15, the range of 5^,2.5 for case 2 is 
well bracketed by that for case 1 between L u x> = 10 and 
100 (in units of 10 51 erg s _1 ). For this reason, we take 
case 1 as our fiducial model (and as an optimum case for 
strong r-processing) in the following sections. 

4.2. Neutron-to-seed Ratios 

The neutron-to-seed ratios at the onset of the r-process 
(dehned at T 9 = 2.5), (Yn/lhKs, are shown in Figure 16. 
We find sizable (Y n /Y h ) 2 . 5 > 1000 (up to - 10 8 ; not 
displayed in Fig. 16) in the innermost region. Except for 
the inner winds, however, (Y n /Yh)2.5 is nearly constant 
with D for L Vt a > 100 (in units of 10 51 erg s _1 ), being 
(Yn/Yh)2.5 ~ 120 at most. This is large enough to expect 
the formation of the third r-process peak (A = 195) and 
beyond, but not fission cycling. A high neutron excess 
(^0,2.5 ^ 0.2) is needed for a strong r-process at the 
modest S (< 100/cb; Figure 7) and r (> 10 ms; Figure 8) 
for the winds except for the innermost region. 
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Fig. 18. — Average atomic mass numbers of heavy nuclei (Z > 2) 
for case 1 (solid lines) and case 2 (dotted lines) as functions of D for 
L, y ,o (in units of 10 51 erg s — 1 ) denoted in the legend. The vertical 
dashed line indicates D = 5i?s. 

4.3. a-Particles versus Heavy Nuclei 

Figure 17 shows the mass fractions of a-particles in 
the final nucleosynthetic abundances (solid lines). In 
the inner winds, the a concentration is significant ow- 
ing to the high S and short r. The r-process thus starts 
from the seed abundances formed i n neutron-rich quas i 
nuclear equilibrium (QSE; A - QO. IMever etaII H998b). 
Note that the a-abundances do not change after charged- 
particle freezeout (Tg ~ 4; Fig. 14). In the innermost 
winds, the extremely high S and short r lead to neutron- 
rich freezeout. This is the reason why Y c _2.5 decreases as 
D approaches the inner tip of the torus (Fig. 15). The 
low mass fraction of a-particles in the outer winds in- 
dicates that the r-process starts from seed abundances 
formed in conditions close to neutron-rich N SE (A ~ 80, 
lHartmann et alj|1985t IWanaio et al.ll2011aft . The nucle- 
osynthetic abundances are thus dominated by heavy el- 
ements, not by a-particles, in the outer winds. 

Figure 18 shows the atomic mass numbers of the fi- 
nal products, Ah, mass-averaged over heavy nuclei with 
Z > 2. We find that A^ is a flat function of D in the 
outer winds for all I^o, similar to (1^/^)2.5 in Fig. 16. 
This suggests that in the outer winds the torus- Y c (as- 
sumed to be constant over the entire range of D) pre- 
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Fig. 19. — Nucleosynthetic abundances for case 1 (top) and case 2 
(bottom), mass-integrated between D = 2R§ and 10-Rs, for the 
values of L^o ( m units of 10 51 erg s — 1 ) denoted in the legend. 

dominantly determines A^, rather than the modest S or 
r, which both exhibit gradients with distance D (Figs. 7 
and 8). In the innermost winds (D < 3i?s), however, the 
high S and short r play crucial roles. Except for the in- 
nermost winds, Ah ranges from 60 to 220, encompassing 
nuclei from the trans-iron to the actinide region, but well 
below the ne utron-induced fissioning p oint (A ~ 290, 
e.g., Fig. 3 in lGorielv fc Clerbauxl[T999l) . This is a con- 
sequence of Y e ,2.5 > 0.17 in the outer winds (Fig. 15), 
which is still too high to expect fission cycling at the 
modest values of S and r. 

4.4. Mass-integrated Abundances 

In order to evaluate the net abundances for each £„.o, 
the nucleosynthetic yields are mass-integrated over the 
entire torus range between D — 2R§ and 10i?s- For the 
j-th L u q, the abundance of nuclide A is calculated with 
Eqs. (13) and (14) as 



(17) 



Figure 19 shows the mass-integrated nucleosynthetic 
abundances for selected L V fl cases. As noted in § 4.3, the 
effect of neutron-induced fission is expected to be negligi- 
ble. In order to roughly include the effect of spontaneous 
and /3-delayed fissions, we simply added the abundances 
with A > 256 (all expected to decay by fission) such as 



o solar r-process abundance z 
case 1 all D < 5 RO < 
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Fig. 20. — Time-integrated nucleosynthetic abundances for the 
entire torus for case 1 (top panel) and case 2 (bottom panel). The 
calculated abundances for case 1 are in good agreement with the 
solar r-process distribution (circles; scaled to match the 3rd and 
2nd peak heights for case 1 and case 2, respectively). The contri- 
butions from the inner (D < 5Rs; green lines) as well as innermost 
(D < 3-Rs; light-blue lines) winds are sub-dominant for case 1, but 
important for case 2. 

Actual abundances will depend on the (highly uncertain) 
decay chains and the abundance distribution of fission 
fragments. The sharp abundance peak at A ~ 140 for 
L„ :0 = 2 x 10 53 erg s -1 (case 1; solid red line in Fig. 19) 
is formed by fission fragments. The effect of fission for 
the other cases are however unimportant. 

The time-integrated yield of nuclei of atomic mass 
number A for the entire torus region is calculated as 



Y A = J^(Y A ) itj Ami, 



(19) 



Y 



A/2 



2Y A 



Y A/2 (A> 256). 



(18) 



making use of Eqs. (15) and (16). In Fig. 20, the re- 
sulting yields are compared with the solar system re- 
process abundances (circles; vertically shifted to match 
the height of the third and second r-process peaks for 
cases 1 and 2, respectively). For case 1, we find good 
agreement of the calculated abundances with the solar 
r-process distribution over the entire range of the re- 
process, A = 90 — 210. No overproduction of the N = 50 
(A ss 90) nuclei can be seen, which were problematic in 
the wind from PNSs born in SNe (|Woosley et al.lll99l 
iWanaio et aT1[200l . This can be understood from the 
fact that the r-processing in the outer winds of BH-tori 
starts from NSE (or a-deficient QSE) seeds (A ~ 80), 
rather than from (a-rich) QSE seeds (A ~ 90). For 
case 2, however, the production of heavy r-process el- 
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ements beyond A — 130 drops steeply. 

Figure 20 also shows the time-integrated abundances 
from the inner (D < 5i?s) and innermost (D < 3i?s) 
torus. For case 1, we find similar abundance curves for 
both regions, but with sizably smaller amounts than the 
total production (10 and 100 times smaller, respectively). 
This indicates that the low mass ejection rates (Fig. 5) 
from the inner and innermost torus diminish the role 
of their nucleosynthetic contributions, regardless of their 
high S and short r. For case 2, on the other hand, the 
contributions from inner and innermost regions dominate 
the abundances heavier than A ~ 140. 

Note that, for case 1, the deficient Pb (A = 206-208) 
and Bi (A = 209) abundances relative to the other r- 
elements for D < 5i?s (Fig. 20) are due to a fundamen- 
tal difference in r-processing between the inner and outer 
winds. The modest r of an outer wind trajectory leads 
to an r-process at high temperature (Tg ~ 1), in which 
case the nucleosynthetic flow approximately follows the 
path determined by the (n,7)-(7,n) equilibrium. In 
the inner region, however, a shorter r results in an re- 
process at substantially lower temperature (Tg < 0.5), in 
which the flow is determined by the competition between 
(n,7) reactions and /3-decays ("cold r-process", Wanajo 
2007; also Blak e & Schramm 1976; Panov & Janka 2009). 
iWanaiol (|2007l ) shows that this non-equilibrium process 
leads to underabundant Pb. 

The ejected mass of r-processed nuclei, defined as those 
with A > 100, is derived as 



m r e j = m 



A>WQ 



'1.30 x 10- 3 M Q (case l).„ n s 
x lO- 3 M (case 2f, ) 



1.16 



: m 



showing a dominance of r-processed material (66% and 
59% for cases 1 and 2, respectively) in the total ejecta 
(Eq. (16)). This is a consequence of the fact that the r- 
process starts from the NSE (or a-deficient QSE) seeds 
in most cases (§ 4.3). The mass of fission products in the 
ejecta can be calculated as well, 

v _ /4.25 x 10- 5 M Q (case 1} 

a 2^ AYA - \ 1.19 x lO- 7 M (case 2},> 

accounting for only 3% and 0.01% of the total r-process 
material for cases 1 and 2, respectively. The mass of 
Eu in the ejecta, the element taken as representative of 
r-elements in Galactic ch emical evolution studies (e.g., 
Ilshimaru fc Wanajolll999T ). is estimated to be 

= J 6.77 x 10- 6 M o (case 



I 3.98 



x 1O- 7 M (. 



v case 



.4=151,153 

Figure 21 shows the ejecta-mass histogram as a func- 
tion of Y Cy 2.5 for case 1 (solid lines) and case 2 (dot- 
ted lines). The color coding in Fig. 21 discriminates all 
abundances, those with all yields having A > 100 (all 
r-abundances) , and those with all yields having A > 180 
(the third peak abundances and heavier). We find that, 
for case 1, all 1^,2.5 values contribute with similar weights 
to the total abundances. For case 2, the values of 
^0,2.5 ~ 0.3-0.4 dominate the total ejecta. Wind tra- 
jectories with low Y e are, however, crucial for the pro- 
duction of the r-nuclei with A > 100 (F C! 2.5 < 0.35), 
and in particular with A > 180 (F C! 2.5 < 0.20). This 
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Fig. 21. — Histogram of the ejecta-mass distribution Am c j (nor- 
malized by the total ejecta mass, m c j) for case 1 (solid lines) and 
for case 2 (dotted lines) as functions of y e ,2.5- Different colors cor- 
respond to all ejecta and those with all yields having A > 100 and 
those with all yields having A > 180, respectively (slightly shifted 
in the horizontal direction for visibility). 

is due to the presence of modest S and r in the ejecta 
(Figs. 12 and 13), which demands relatively low Y e val- 
ues for strong r-processing. Figure 21 also shows that the 
later wind outflow with L v < 10 52 erg s _1 (t > 100 ms 
and Y c ~ 0.5), not included in our calculations, would 
make no relevant contribution to the r-abundances. 



5. BH-TORUS WINDS AS THE ORIGIN OF r-ELEMENTS 

In § 4, we found that our fiducial model (case 1) of 
BH-torus winds leads to the full r-process with the solar- 
like r-pattern for A = 90-210. The ejecta mass of the 
r-processed matter (A > 100), m r , c] = 1.30 x 1O~ 3 M 
(Eq. (20); case 1), is more than a factor of 10 larger than 
that needed for CCSNe to be the dominant source of 
the r -process elements in the Galaxy (Mathew s &: Cowanl 
1990). More specifically, the ejected mass of Eu (a 
nearly pure r-process element), TOEu.oj — 

6.77 x 1O~ 6 M 

(Eq. (22); case 1), is a factor of 60 greater than that 
needed for CCSNe to be the majo r source of the Galac- 
tic Eu (|Wanaio &: Ishimaru 200H). Provided that our 
fiducial model (case 1) with Mbh = 3M and a torus 
mass of Mtorus ~ O.1M (Eq. (11)) is representative of 
NS-NS (or BH-NS) mergers, the canonical CCSN event 
rate of ~ 10~ 2 yr _1 implies that a time-averaged Galac- 
tic merger rate of ~ 2 x 10~ 4 yr _1 would be needed. This 
is consistent with the upper bound of a merger event rate 
between 7 x 10~ 6 and 3 x 10~ 4 yr _1 derived by popu- 
lation synthesis methods for NS- NS and BH-NS b i naries 
(at solar-metallicity conditions, Be lczvnski et al.l I2002D 
and of statistical result s based on binary pu lsar surveys 
(2 x 10~ 5 -3x 10~ 4 yr- MKalogera et"aTll2004f ). For case 2 
(with inclusion of reactions of Eqs. (l)-(4) at Tg > 10), 
we find inefficient production of r-elements beyond A ~ 
130. However, a similar constraint to the Galactic merger 
rate may be applied for this case, owing to the simi- 
lar r-processed ejecta mass of m rci = 1.16 x 1O _3 M 
(Eq. (20); case 2). 

It is important to note that the actual ejecta mass 
of r-abundances could be higher than that we obtained 
on the basis of our simplified BH-torus wind model. 
The total ejecta mass, m c] w 2 x lO~ 3 Af (Eq. (16)), 
can be even 10 times smaller than that expected from 
hydrodynamical simulations of compact object mergers 
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Fig. 22. — Time-integrated nucleosynthetic abundances for the 
entire torus for case 1 (solid line) and case 2 (dashed line). Cases 1 
and 2 arc compared with the spect roscopic abun dances of r- 
enhanced (CS 31082-001; filled circles, Sncdcn ct al. 2003) and r- 
deficient (HD 122563; open circles, Honda ct al. 2006) Galactic 
halo stars, respective l y. Fo r HD 122563, the Cd and Lu values are 
from Rocdcrcr ct al. (2010) and the Ge value is from Cowan ct al. 
(2005). For both stars, the abundances are vertically shifted to 
match the calculated Eu abundances. 

(but not of the subsequ ent BH-torus accre tion with its 
neutrino-driven outflow; Uanka et al.l 119991 ). There will 
be additional contributions from the early mass loss 
due to the tidal ejection of neutron-rich matter during 
the m erger event (|Freiburghaus et aT1ll999tlGorielv et al.l 
2005) and neutri no-driven outflows from HMNSs (for 
the NS-NS case, D essart et al.l I2009D or /and magneti- 
cally driven outflow s from HMNSs (jRezzolla et aLll2011l: 
Shiba ta~et al.ll201lj) . Moreover, the centrifugal force due 
to rapid rotation, which we do not take into account, ef- 
fectively reduces the gravity from the central BH. This 
could lead to la rger mass ejection rates than those ob- 
taincd here (see lOtsuki et al.ll2000t IWanajo et aT1l200lL 
for PNS winds). Magnetic fields, which are not consid- 
ered in our models either, could play a significant role for 
the mass ejection, producing viscously driven or MHD- 
driven outflow s in addition to the considered neutrino- 
driven ejecta (M etzger et al.l 120091 ) . It should also be 
noted that the production of heavy r-process elements 
beyond A ~ 130 is marginal if we consider our case 2 to 
be representative of NS-NS and BH-NS mergers. Case 2 
takes into account v Cl D c , e~ , and e + captures on free nu- 
cleons also for Tg > 10, resulting in significantly higher 
minimal 1^.2.5 (~ 0.3) than that in case 1 (~ 0.2). How- 
ever, the detailed evolution of Y c will be highly depen- 
dent on the density structure and the uncertain neutrino 
field in the vicinity of the BH-torus, while we consider 
a simplified spherical wind model and do not have neu- 
trino transport results for BH-tori at bound. Effects of 
magnetic fields noted above would also modify the evolu- 
tion of Y c . General relativistic effects (which we consider 
only in the framework of spherically symmetric configu- 
rations) can a lso be important for d etermining the exact 
Y c evolution dCaballero et al.l 120111) . For all these rea- 
sons, the nucleosynthetic outcome of this study should 
be regarded only as suggestive. Future detailed multi- 
dimensional simulations of NS-NS (or BH-NS) merging, 
including the later BH-torus wind phase, are required for 
more quantitative results. 

Nevertheless, the good agreement of our nucleosyn- 
thesis result for case 1 with the solar r-pattern is en- 



couraging. Such an agreement is also important when 
one considers compact binary mergers as the origin of 
early Galactic r-elements in metal -deficient stars wit h 
their uniform (solar- like) r-patterns (jSneden et al.ll2008h . 
In Figure 22, the nucleosynthetic abundances for case 1 
(solid line; as a function of atomic number) are com- 
pared with an r-process enhanced star CS 22892-052 with 
[Fe/H] « — 3.1 (filled circles; scaled to m atch the calcu- 
lated Eu abundance, iSneden et al.ll2003D . We find quite 
a good agreement of our result with the stellar abun- 
dances distribution of CS 22892-052. The good agree- 
ment of model and observed abundances is obtained if 
our BH-torus wind model represents typical NS-NS (or 
BH-NS) merger events, if the minimal 1^,2.5 (at the on- 
set of r-processing) is as low as in case 1 (< 0.2), and if 
the evolution of L Vb and L Pe evolutions in Eqs. (5) and 
(6) is appropriate (e.g., a in Eq. (6) is not very differ- 
ent from unity and /3 = 1.3 and the factor of 2 within 
the bracket of Eq. (6) are good choices). These should 
be tested by future long-time hydrodynamical simula- 
tions of NS-NS and BH-NS merging events and of the 
evolution of their remnants. Note that, even if the pro- 
duction of r-elements heavier than A ~ 130 is marginal 
as found in case 2, our model could be a possible ex- 
planation for the abundance distribution found in an r- 
process deficient star HD 122563 with [Fe/H] w -2.7 
(open circles; scaled to m atch the calculated Eu abun- 
dance, iHonda et all 120061 ). This suggests that NS-NS 
and BH-NS mergers could at least be the origin of some 
trans-iron elements up to Z ~ 50 (^4 ~ 120). 

For an observational support, future modeling of the 
Galactic r-element evolution that is not in conflict with 
spectroscopic results of meta l-deficient stars will be im- 
portant (e.g. iPrantzosI 120061 ). In addition, radioactive 
decays of r -processed ejecta can lead to faint SN-like 
transients (iLi fe Paczvhskil 119981 : iMetzger et all l201ti 
iRoberts et al.l 1201 lh . Future detections of such signals, 
possibly as accompanying events of short GRBs, will 
be a direct support for the occurrence of r-processing 
in NS-NS and BH-NS m e rgers. According to the es- 
timates bv IMetzger et all ((2010), a peak luminosity of 
~3x 10 41 erg s" 1 could be expected from the BH-wind 
ejecta with m r , c - s ~lx 1O _3 M (Roberts et al. (2011), 
however, obtained only slightly higher R-band luminosi- 
ties for considerably higher ejecta masses). It will be dif- 
ficult to distinguish the wind component from the early, 
tidally ejected counterpart. Both could contain similar 
amounts of radioactive matter. However, the velocities 
as well as the geometry of the r-processed ejecta might 
be distinctive between them. Our result implies that 
the r-processed BH-wind ejecta reach out widely from 
the torus with velocities of ~ 0.1c (except for the inner- 
most region, Figs. 6 and 11). In contrast, tidal ejecta 
come from the tips of (one or two) spiral arms, and are 
expected to achieve larger outgoing velocities (> 0.2c, 
IRoberts et alJl20ilI ). These could lead to some different 
features in the light curves (e.g., line profiles or the time 
of the peak), potentially distinguishable by observations. 

6. SUMMARY 

We investigated r-process nucleosynthesis in the 
neutrino-driven outflows from a BH accretion torus 
(BH-torus winds) formed in NS-N S (or BH-NS) merg- 
ers. Different from previous works (|Surman "eTaLl 120081: 
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iMetzger et aL|[20loh we investigated a time- and space- 
dependent semi- analytic model setup. The BH-torus 
wind models were constructed by considering spherically 
symmet ric, general relativis tic neutrino-driven wind so- 
lutions (jWanajo et al.1 120011 ) with the gravitational po- 
tential of a central BH of M B h = 3M . The BH ac- 
cretion torus (around a rotating BH) was assumed to 
lie between D = 2R S and 10i? s (-Rs = 8.86 km) from 
the center. Each wind trajectory reaching away from 
the torus was obtained by assuming a hypothetical neu- 
trinosphere in the spherical wind model with the radius 
R„ = D. 

In the innermost wind region (D ~ 2i?g), the effi- 
cient energy deposition due to vv pair annihilation to 
e~e + pairs leads to very high entropies (S > IOO-IOOO/cb 
per nucleon) and short expansion timescales (r < 1— 
10 ms). This allows for a strong r-proce ss regardless 
of Y c (even with Y c > 0.5 when r < 1 ms, IMeverll2002j) 
or no production of heavy elements (in the relativistic 
wind s with S > 1000fc B , iLemoind l200l iBeloborodovl 
2003). However, the small mass ejection from the in- 
nermost torus makes this contribution to the total nu- 
cleosynthetic abundances essentially negligible. In the 
outer wind regions (D > 5i?s), on the other hand, 
the dominant heating is due to v e and v e captures on 
free nucleons as in the case of PNS winds, resulting in 
modest entropies (S ~ 30fce) and expansion timescales 
(t ~ 100 ms). Low Y c values (< 0.2 at T 9 = 2.5) are thus 
essential for a strong r-process. The contribution from 
outer wind regions dominates the total nucleosynthctic 
abundances because of their greater mass ejection rates. 
Note th at our BH-torus model does not pred i ct a vp- 
process (jFrohlich et all [20061: IPruet et all 120061: iWanaiol 
I2006a|) as s uggested in the case of (proton-rich) collapsar 
disk winds (jKizivat et al.ll2010l) . 

The mass-integrated nucleosynthetic abundances, ob- 
tained with a phenomenological time evolution of neu- 
trino luminosities, are in good agreement with the solar 
r-pattern over the entire r-process range of A = 90 — 210, 
when the neutrino-matter interactions are considered 
only for Tg < 10 (case 1). This can be taken as the op- 
timal case for strong r-processing in our models. How- 
ever, when 14, 24, e~, and e + captures on free nucle- 
ons are taken into account also for Tg > 10 (case 2), 
the production of heavy r-elements beyond A ~ 130 



drops off steeply. The total ejecta mass of the BH- 
torus wind of our simplified model is calculated to be 



m j » 2 x 10 3 Mq, in which the r- processed matter 



dominates (m r , e j > 1 X 10 M©). Provided that our 
BH-torus wind model is representative, NS-NS and BH- 
NS mergers can produce all (or at least a significant part 
of) the Galactic r-abundances if the event rate (aver- 
aged over the Galactic history) was ~ 2 x 10~ 4 yr" 1 , 
which is consistent with the upper bound of present 
(solar-metallicit y) population synthesis results (7x 10~ 6 - 
3 x 10~ 4 yr-\ IBelczvnski et all I2002D and of statisti- 
cal results based on binary pulsar surveys (2 x 10 -5 - 
3 x 10" 4 vr" 1 . Kalogera et al.ll200l . This implies that 
BH-torus winds from NS-NS and/or BH-NS mergers 
could be major (or partial) production sites of the r- 
process elements in the Galaxy. It should be noted that 
the actual ejecta mass of r-abundances could be sub- 
stantially higher than the estimate based on our BH- 
torus wind model, in which we do not consider any other 
effects than neutrino-driven outflows, namely, centrifu- 
gal force, magnetic fields, or tidally ejected neutron-rich 
matter from NS disruption. 

Obviously, more elaborate hydrodynamical studies of 
the considered astrophysical site are needed to obtain 
information of the neutrino field that controls the dy- 
namics as well as the neutron-richness in the BH-torus 
winds. Note that NS-NS and BH-NS mergers are also 
suggested to be the sources of short GRBs. An in- 
teresting possibility in this context is that the radioac- 
tive neutron-rich nuclei ejected during and after mergers 
migh t lead to detectable transient electromagnetic sig- 
nals dLi fc Paczvhskilll998l; Kulka rni 2005; IMetzger etall 
I20l5 iRoberts et al.|[201l|) . Studies of Galactic chemical 
evolution will also be important to test the contributions 
of NS-NS and BH-NS mergers to the enrichment history 
of the r-process elements, in particular their role in the 
early Galaxy. 
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